Spatial correlation of solar wind fluctuations and solar cycle dependence 
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We investigate the spatial correlation properties of the solar wind using simultaneous observations 
by the ACE and WIND spacecraft. We use mutual information as a nonlinear measure of correlation 
and compare this to linear correlation. We find that the correlation lengthscales A of fluctuations in 
density p and magnetic field magnitude \B\ vary strongly with the solar cycle, whereas correlation 
lengths of fluctuations in B field components do not. We find A(|i3|) ~ f 20-Rb and A(|-B|) ~ 270Re 
at solar minimum and maximum respectively and A(p) ~ 75Re and A(p) ~ f 70-Rb at minimum and 
maximum. The components of the B field have A -^ A(|B|) at minimum. 

PACS numbers: 



I. INTRODUCTION 



In situ solar wind plasma observations over the last 30 
years show that its local properties at 1 a.u. are modu- 
lated by the solar cycle [1] and that this has an impact on 
the magnetosphere and ionosphere Q ■ There is evidence 
for more mixing of fast and slow solar wind plasma at so- 
lar maximum [3| and changing magnetic for complexity 
over the cycle [4 1 Q • One method for quantifying the ef- 
fect of solar activity on the solar wind is to determine the 
correlation length scale of fluctuations measured in-situ. 
In the past this was done using the Taylor hypothesis [q| 
and long timeseries from single spacecraft (e.g. Bruno et 
aZ Q ) • More recently, multiple spacecraft studies have 
used linear correlation to calculate a typical correlation 
length for magnetic field in the solar wind. The ACE and 
WIND spacecraft give an excellent range of separations 
and enable long baseline correlations to be calculated for 
magnetic field vector B and magnitude \B\ as well as 
solar wind bulk velocity and density p. These observa- 
tions have been used to show anisotropy in the solar wind 
and measure the Taylor microscale Q- Quantitative 
knowledge of the spatial correlation of the fluctuations of 
the flow is relevant to attempts to understand the solar 
wind in terms of locally evolving turbulence [10| 11 ill . It 
also has implications for cosmic ray propagation |12l | [l3| 
which shows solar cycle dependence jl4| . 

Previous studies used linear correlation measures and 
did not investigate solar cycle effects on the measure- 
ments of correlation length. Here we use a nonlinear 
measure of correlation, mutual information [^ [l5| [l6[ 
[17| p^ . alongside linear correlation coefficient. We cal- 
culate correlation lengths for the components of B, as 
well as |i?| and p. We also investigate the behaviour of 
the correlation length of these quantities over the recent 
solar cycle, using ACE and WIND spacecraft data. 



II. MEASURING CORRELATION 

We use two methods for estimating correlation: corre- 
lation coefficient as a linear measure, and normalised mu- 
tual information as a nonlinear measure. The linear cross 
covariance provides a measure of correlation between two 
signals A and B, defined by: 



C{A,B) 
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(1) 



where E[. . .] denotes the mathematical expectation value 
and A = E[A]. Mutual information (MI) quantifies the 
information content shared by two signals A and B. For 
discrete signals we can write the MI as: 



/(A,B)=^P(a„6,)log2 



P{ai)P{b/) 



(2) 



Here the signals A and B have been partitioned into 
an alphabet (an exhaustive discrete set which spans 
the possible values the signal can take) so that A — 
{oi, . . . jQi, . . .am} where oi and am are the extrema of 
A found in all data considered. The discretised signal 
takes value a^ with probability P{ai) and similarly for bi 
we have P{bi), while P(ai,bj) is the joint probability of 
ai and bj . The chosen base of the logarithm defines the 
units in which the MI is measured. Normally base 2 is 
used, so that the MI is measured in bits. If one defines 
the entropy of a signal as 



H{A) - -^P(a.)log2(P(a,)), 



(3) 



then MI can be written as a combination of entropies [l| 
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I{A,B) = H{A) + H{B)-H{A,B) 



(4) 



The calculation of the entropies needed to form the 
MI is not trivial, as there is some freedom in the method 



of discretisation of the signals and in the method used 
to estimate the probabihties P{ai), P{hj) and P{ai,bj). 
There are many different methods currently used, sum- 
marised and compared by Cellucci et al. 19] and Kraskov 
et al. [20]. For spacecraft observations of the solar wind, 
we use a discretisation based on the standard deviation a 
of the data. Only data within 5cr of the mean is consid- 
ered and bins of ^cr are used, giving 20 bins in total. The 
form of MI used here is the normalised mutual informa- 
tion (NMI) [U, which is I{A, B) normalised by the joint 
entropy H{A, B) so as to remove the dependence on the 
entropy of the solar wind at the time of observation: 
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This gives results in the range < NMI < 1 facilitating 
direct comparison between different periods of data. 



III. RESULTS 

A. The Datasets 

To investigate the effect of the solar cycle on spatial 
correlation in the solar wind we take data from periods 
as near to solar activity minimum as possible, that is in 
1998, 2005 and 2006, and from solar maximum in 2000. 
In 1998 WIND was returning towards Earth from the 
Sun-Earth libration point while ACE orbited the libra- 
tion point; in 2005 and 2006 both ACE and WIND or- 
bit the libration point. 17 periods of data were chosen 
with different spacecraft separations, giving a total of 48 
days of data. At solar maximum in 2000 WIND is or- 
biting the Earth and therefore only short periods, when 
the spacecraft is suitably far from the magnetopause and 
bow shock, can be used. Four such periods were used giv- 
ing a total of 15 days of data. Two minute cadence data 
are used for all variables. We use two different lengths 
T for our measurement window to optimise for the two 
distinct timescales in the power spectrum. We choose 
tl — 960 minutes to investigate the large scale struc- 
tures and Ts = 200 minutes to access the inertial range 
of turbulence in the solar wind. A shorter window can- 
not be used as the data in the window then becomes 
too short for a reliable correlation estimate to be made. 
NMI ([5]) and correlation coefficient ([T]) are calculated for 
windows rg and tl as they are moved along the data. A 
value for NMI and correlation coefficient is thus obtained 
every two minutes for each window length, 24 hour av- 
erages are calculated and plotted against separation in 
Figures [1] to [H 



B. Analysis 

Figure [T] plots linear correlation and NMI versus space- 
craft separation for p , \B\ and B components at solar 
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TABLE I: Values of correlation length A using window tl, 
with 95% confidence bounds, in Earth radii, calculated by 
least squares fitting of exponentials to the 24 hour mean NMI 
and linear correlation coefficient during solar minimum (1998, 
2005, 2006) and solar maximum (2000). Those values marked 
with * are fits that resulted in an R^ value of less than 0.1. 



minimum using r^. Following [8| we fit an exponential 
function y = aexp(a;/A) to the data using nonlinear least 
squares, with a = 1 for correlation coefficient. The val- 
ues of A calculated are shown in Table HI The quality of 
the fit is indicated by the error on the correlation length, 
calculated as the 95% confidence bound of the nonlinear 
least squares fit. The errors on these values arise pri- 
marily from the scatter of the measurements due to the 
highly variable nature of the solar wind. 
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FIG. 1: NMI and linear correlation for tl = 960 minutes, 
24 hour averages plotted here versus spacecraft separation. 
Density p {{+) and dashed line) is compared to magnetic field 
magnitude |S| ((•) and solid line) in the two left hand panels. 
In the right hand panels the y and the z components (GSE) 
of B are compared. The lines are inverse exponential fits to 
the data. 

We anticipate the spatial variation of correlation coef- 
ficient to be of exponential form. However NMI, which is 
a nonlinear measure, may have more complex behaviour; 
an exponential function provides a simple parameterisa- 
tion of the data. The NMI gives p a longer correlation 
length than \B\, whereas the correlation coefficient gives 
the reverse. Both give a lower correlation for p than \B\. 
If both p and \B\ NMI fits are set to have the same in- 
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FIG. 2: 24 hour mean NMI and correlation coefficient coeffi- 
cients, between ACE and WIND \B\ and p for window length 
tl = 960 minutes, for selected periods in 1998, 2005 and 2006 
(solar minimum: (•) and solid line) and 2000 (solar maxi- 
mum: (-I-) and dashed line). The lines plotted are inverse 
exponential fits to the data. 



tercept then the NMI gives similar results to those found 
with the cross corelation, however as exact form of the 
NMI dependence on separation is unknown these corre- 
lation lengths remain a local estimate. 

TableUalso shows that at minimum the y and z compo- 
nents of B have similar correlation lengths to each other 
and |i3|. The scatter on the correlation calculated for the 
components of B and |i3| is larger than for p. The corre- 
lation coefficient gives a clearer decay and the estimates 
of correlation length are correspondingly more accurate. 

We now consider variation with solar cycle. Figure [2] 
plots linear correlation coefficient and NMI against space- 
craft separation for \B\ and p. It shows that both the 
linear and the nonlinear measure of correlation yield a 
higher spatial correlation in both p and |i?| at solar max- 
imum. To quantify this we again fit an exponential of 
the form y — aexp{x/X), with a — 1. This provides 
a reasonable fit in all cases except for the NMI at so- 
lar maximum for p, which also shows the largest scatter. 
The calculated values of correlation length A are shown 
in Table HI We see that the correlation length measured 
at solar minimum is systematically smaller than that at 
maximum. The value at maximum for |-B| is within error 
of previous estimates for the same interval Q . 

Figure [3] shows the spatial correlation of components 
of B. The points from solar maximum and minimum do 
not show distinct behaviour, unlike p and \B\ in Figure 
[2l Fitting an exponential function to the data, as before, 
yields the values of the correlation lengths for the com- 
ponents of B shown in Table HI The By and B^ values 
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FIG. 3: 24 hour mean NMI and correlation coefficient coef- 
ficients, between ACE and WIND By and Bz using window 
length Tl ~ 960 minutes, for selected periods in 1998, 2005 
and 2006 (solar minimum: (•) and solid line) and 2000 (so- 
lar maximum: (-I-) and dashed line). The lines plotted are 
inverse exponential fits to the data 
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FIG. 4: 24 hour mean correlation coefficient coefficients, cal- 
culated between ACE and WIND |_B| and p measurements 
at solar minimum in 1998, 2005 and 2006 (top panels) and 
solar maximum in 2000 (bottom panels). The (•) and solid 
line fits correspond to tl = 960 min and the (-I-) and dashed 
line fits correspond to ts — 200 min. The lines are inverse 
exponential fits to the data. 



distinct from the values for \B\ and p. 

Finally, Figure 3] gives a comparison between the mea- 
sured correlation coefficient from the two different time 
windows, tl = 960 minutes and ts = 200 minutes. The 
shorter window is within the inertial range of solar wind 
turbulence; it gives slightly smaller values for the corre- 
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lation length (A™„(|-B|) 

Re) than the longer time window {Xmini\B\) — 118 _i5, 

Amm(p) = 75 -g). This difference is more obvious at 
solar maximum where the correlation lengths measured 
by the smaller window are Amaa:(|-B|) 
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for \B\ and Xmaxip) = 116 -22 Re for p; compared to 
Xmaxi\B\) = 274 t^'^ Re and \,nax{p) - 167 ±I§ Re 
using the longer window tl- Large coherent structures 
within the solar wind may contribute to the change in 
correlation length seen on timescales tl above the iner- 
tial range value Tg; however the difference is not large. 



IV. CONCLUSIONS 

We have used simultaneous data from the ACE and 
WIND spacecraft separated by between 30 and 220Re, to 
calculate the spatial correlation of solar wind p, |-B| and 
B components. Cross correlation and normalised mutual 
information were used as alternative linear and nonlinear 
measures of correlation. At solar minimum we have 48 
days of contemporaneous WIND and ACE observations. 
We determine the correlation lengthscalc A(p) for p to be 
A(p) = 75Re. This is smaller than X{\B\) by - 1.6. For 
the components of B, we find X{Bz) is within error of 
A(|i?|) and X{By) is slightly above the range of errors. 

At solar maximum we have 15 days of contemporane- 
ous WIND and ACE data. This allows us to investigate 



the effect of the solar cycle. We find that A(|-B|) and A(p) 
at solar maximum are larger than at minimum by ~ a 
factor of 2. In contrast the components of the magnetic 
field show weak variation with solar cycle. 

We have used two window lengths; one within the in- 
ertial range (200 min) and the other on longer timescales 
(960 min). The window on inertial range timescales gives 
values of A, that are systematically shorter, but within 
the errors, of the values from the larger window. The 
solar cycle dependence of the correlation length is in all 
cases independent of our chosen window sizes. 

Our result that p and \B\ show variation with solar 
cycle, whereas the components of B do not, is consistent 
with the idea that the correlation in \B\ and p is, at least 
in part, of solar origin. Intriguingly this behaviour per- 
sists when we restrict our analysis to timescales within 
the inertial range. This is consistent with recent single- 
spacecraft results Q which show solar cycle depen- 
dence within the inertial range of magnetic energy den- 
sity fluctuations. 
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